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ABSTRACT 

We propose an estimator of the microlensing optical depth from pixel 
lensing data that involves only measurable quantities. In comparison to the 
only previously proposed estimator, it has the advantage of not being limited 
to events with large magnification at maximum, and it applies equally well to 
satellite and ground-based observations. 



Subject headings: gravitational lensing 
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1. INTRODUCTION 

An important physical quantity that we want to determine from gravitational 
microlensing data is the microlensing optical depth. This is the number of gravitational 
lenses lying within one Einstein radius around the line of sight to a given source. 

When monitoring individual stars as proposed by Paczyhski (1986), the optical depth is 
estimated from observational data as a weighted sum of the Einstein times of the individual 
events (see eq. J7J below; Griest 1991). The Einstein times are determined by fitting the 
theoretical microlensing lightcurve to the measured lightcurve. The fit is almost degenerate 
in the three fit parameters Einstein time, maximum magnification, and star flux in absence 
of lensing (Gould 1996; Wozniak & Paczyhski 1997). In particular, the fitted value of 
the Einstein time depends sensitively on the value of the star flux in absence of lensing, 
especially at large maximum magnifications. 

In pixel microlensing, which is gravitational microlensing of unresolved stars, the star 
flux in absence of lensing is not generally accessible, since the combined light of all the stars 
contributing to a pixel is monitored simultaneously. Thus, because of the above-mentioned 
degeneracy, the Einstein time cannot be reliably determined from a fit to the lightcurve. In 
this case, the optical depth cannot be estimated by means of the usual expression, and this 
seems to hinder the possibility of obtaining it from pixel lensing data. 

A partial solution to this problem has been proposed by Gould (1995 & 1996) under 
conditions that may be possible in future satellite observations but are definitely not 
typical of present microlensing searches towards the Andromeda galaxy M31 and the 
Large Magellanic Cloud (Crotts 1992, Baillon et al. 1993, Crotts and Tomaney 1996, 
Tomaney and Crotts 1996, Ansari et al. 1997 for M31; Melchior et al. 1997 & 1998 for 
the LMC). He assumed regular observations with constant exposure times, a constant level 
of noise dominated by photon counting, and large magnifications at maximum so that 
the microlensing lightcurve assumes a completely degenerate form (eq. || below). The 
present pixel-lensing searches instead have irregular time gaps between observations, show 
a variable level of noise depending on factors other than photon counting, and expect 
relatively low magnifications for which the microlensing curve is not well approximated by 
the completely degenerate form. 

In this letter, we present an estimator of the optical depth in pixel microlensing 
that applies to all magnifications, involves measurable quantities only, and is of direct 
applicability to the data. It is based on the substitution of the Einstein time with the 
product of the flux increase at maximum and of the full event duration at half-maximum, 
defined as the total time during which the flux increase exceeds half of its maximum value. 
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2. USUAL OPTICAL DEPTH ESTIMATOR 

The flux increase in a microlensing event is given by (for pointlike lens and pointlike 
source; Einstein 1936) 

AF(t) = Ff(u\t)), (1) 

with 

m = - 1. (2) 

Here F is the source flux in absence of lensing, u = 0/9 e , 9 is the angular separation of lens 
and source, and 6 e is the angular Einstein radius 
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where m is the lens mass, and s and / are the distances of the source and the lens from the 
observer. 

The relative motion of lens and source during the event is usually well approximated 
by a uniform motion, for which 

/t-t ^ 2 
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u\t)=F+(^-j^) . (4) 

tmax is the time of maximum magnification, (3 is the impact parameter in units of the 
Einstein radius, and t e is the Einstein time. This we define as the time the lens would take 
to cross the radius of the Einstein ring, t e = 9 e / fx, where ji is the relative proper motion of 
lens and source. 

In the limit of large magnifications at maximum, the microlensing curve assumes the 
completely degenerate form 
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(F,f3 -> 0; F/P = const). (5) 



The three parameters F, (3 and t e have reduced to two: F/P and pt e . A fit to the data can 
determine neither the star flux F nor the Einstein time t e separately from the unknown 
impact parameter P (see e.g. Wozniak & Paczyhski 1997). 

The microlensing optical depth is the number of gravitational lenses lying within one 
Einstein radius 9 e around the line of sight to a given source, 



= f 7i9 2 J 2 n{l)dl, (6) 
Jo 
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where the integral is along the line of sight to the source at distance s from the observer, 
and n(l) is the lens number density at distance I along the same line of sight. 

When monitoring N s sources for a time T, the optical depth is usually estimated as 

T = 2N T ^ fit V ^ 

ziv s J events H 1 ^ 

Here r is the optical depth averaged over the monitored sources, the sum is over the 
observed events, and e(t e ) is the detection efficiency (normalized to impact parameters 
< 6 e ) as a function of Einstein time t e . The Einstein time of each event is determined 
experimentally by a fit of the theoretical microlensing lightcurve to the measured lightcurve. 

This usual estimator is based on the relation 

tJNJT = ?-N s t. (8) 

i9<1 2 

In the right hand side, f = N" 1 / rdN s ; in the left-hand side, dN s is the differential number 
of monitored stars, dT is the event rate per source, dT = (2r/7rt e ) df3, and the integration is 
limited to impact parameters f3 < 1. Eq. fl7|) follows from solving eq. (§) for f and rescaling 
the observed number of events dN ohs to the expected number of events dN = TdN s dT by 
means of the detection efficiency e(i e ) normalized to /3 < 1. Explicitly 

y Jl^ = f J^dN ohs = I t e dN = —TN S T. (9) 



3. NEW OPTICAL DEPTH ESTIMATOR 

We now prove a relation similar to eq. (Jjjfy but not involving the Einstein time t e , which 
is very poorly measurable in pixel lensing. Let tf w h m be the event duration at half-maximum 
and A max be the flux increase at maximum. Both are easily determined from a fit to the 
lightcurve. Then 

J t iwhm A max dN s dT = —IF s Tp. (10) 

Here F s is the total flux of the monitored sources, / is 2.1465, and Tp is the optical depth 
averaged over the monitored sources, each source weighted by its flux, Tp = F^ 1 J FrdN s . 
This flux-weighted optical depth Tp equals the averaged optical depth r whenever the 
source luminosity function is, or can be considered as, independent of source position, 
because in this case the averages over source positions and fluxes are independent. This is 
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the case for a single source population.^] 

The proof of eq. ( TO ) is as follows. We express the half-maximum duration tf wnm as 
a function of t e and f3 by finding the value of the impact parameter f3 w at which the flux 
increase equals half of its maximum value. We also express the maximum flux increase 
A max in terms of the source flux F in absence of lensing and of the impact parameter (3. 
We find t fwhm = t e w(P) and A max = F5((3), with w((3) = 2 - (3 2 ) l/ \ (3 2 w = 2 f (5(13)), 
5((3) = f(/3 2 ), and f(x) given in eq. (0). The integral over (3 in eq. (|10D is then 

w((3)5((3)dp = 2.1465. (11) 



o 

The remaining integral over the sources gives the total source flux F s = J FdN s times the 
flux-weighted optical depth r>. Putting the factors together gives eq. fllOf). 

If all events could be detected, then eq. fllQ|) would immediately give an estimator for 
the optical depth, namely 

r F = 2JTF ^ ^fwhmA max (no noise, 100% efficiency). (12) 

s all events 

Notice that apart from the numerical factor it/ 21 ~ 0.73, Tp equals the sum over events 
of the products of the fractional flux increase at maximum A max /F s and of the fraction 
of observing time spent above half- maximum tufom/T. The estimator in eq. ( |T^) does not 
depend on the source luminosity funciton. 

In practice, there is always a minimum detectable A max , if only because of Poissonian 
fluctuations in photon counting. Moreover, the irregularity and the discreteness of the time 
coverage limit the range of event time scales a microlensing search can access. 



We can restrict the sum in eq. (0) - and the integral in eq. (10) - to events with 



A max > A, a convenient detection threshold. We obtain 



TT 



^ . ^fwhmA max , (13) 



events 
Amax > A 



2ITF cS (A) 
where we have defined an effective flux 

F eff (A) = r l J dF<P(F)F^ A/F d[3w((3)6(P), (14) 



lr rhe flux- weighted optical depth may actually find an interesting application for individually monitored 
sources, for example in the microlensing data towards the galactic bulge: comparing it to the unweighted 
optical depth may give information on a spatial segregation of the star populations. 
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4>(F) being the source luminosity function and (3a/f — 2f(2A/F). This is the basic form of 
the new optical depth estimator that we advocate for pixel lensing. Eq. (|T3| ) is written for a 
single source population; for several populations, similar equations apply to each population 
separately 

In the next section we discuss the dependence of the effective flux F eff (A) and of the 
new optical depth estimator on the source luminosity function. 

4. EFFECTIVE FLUX AND SENSITIVITY 

The effective flux F e g(A) is smaller than the total flux F s and is equal to F s when 
A = 0. It depends in principle on the details of the source luminosity function, but in 
practice the ratio F e ff(A)/F s is an almost universal function of the fluctuation magnitude 
m (defined in Tonry & Schneider 1988). This is exemplified in fig. 1 for four luminosity 
functions: (1) the local neighborhood in the K-band (Jahreifi & Wielen 1998); (2) 47 
Tucanae in the K-band (Hesser et al. 1987); (3) the M31 bulge in the i^-band (Rich, 
Mould, and Graham 1993); (4) the galactic bulge in the J-band (Terndrup, Frogel and 
Whitford 1990; Holtzman et al. 1998). The upper panel shows // e g — fi as a function of 
m t hr — rn, where /i e ff is the surface magnitude corresponding to the effective flux, \i is the 
actual surface magnitude, and m t h r is the magnitude corresponding to the threshold flux 
increase A. The lower panel shows the fractional differences in effective flux with respect 
to the mean over the luminosity functions considered, which is also the relative error in 
the estimated optical depth. At very low thresholds (m t h r 3> rn), the effective surface 
magnitude tends to the actual surface magnitude, and this for all luminosity functions. 
At very high thresholds (m thr <C m tip , the magnitude of the brightest star), the difference 
between effective and actual surface magnitudes is simply /i c g- — /z = m — m thr — 0.52, and 
this again for all luminosity functions. The largest deviations from a universal behavior 
occur at detection thresholds close to the fluctuation magnitude, and for the luminosity 
functions considered here they are never larger than 30%. If the statistical errors in the 
determination of the optical depth will be smaller than these systematic errors (we recall 
that the present error on the optical depth towards the LMC is ~ 50% [Alcock et al. 1997]), 
the systematic errors can of course be reduced by measuring the source luminosity function 
in deep exposures, even independently of the microlensing search itself. 

A given microlensing search is only sensitive to a range of Einstein times and radii, and 
can only measure the contributions to the optical depth coming from this range. Indeed, 
with a finite time coverage, there is a minimum event duration set by the finite exposure 
time and the interval between successive images, and a maximum event duration set by the 
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total time span of the observations. Moreover, for sources with finite angular extensions, 
there is a reduction in flux increase with respect to the pointlike case when the lens transits 
the source disk. As a consequence, the fraction of observable events is different according 
to the Einstein time and radius of each event. 

We write the optical depth observable in a given microlensing search as 

r bs = J ir6ll 2 n(v T ,l)s(t e ,9 e )dv T dl (15) 

(compare with eq. ||, where n(l) = J ti(vt, l)dvr)- Here vt = l/J> is the relative transverse 
velocity of lens and source. 

The function s(t e , 9 e ) represents the sensitivity of a given experiment to the different 
contributions to the optical depth, and depends on the observing conditions and data 
analysis and on the choice of the optical depth estimator. We call s(t e , 9 e ) the optical depth 
sensitivity. For the usual estimator, s(t e , 6 e ) equals the detection efficiency as a function of 
t e and 9 e . For the new estimator in eq. ([Hp, s(t e ,9 e ) is 

s(t e , 9 e ) = } ... / F<f>{F) w(/3)6{/3) e obs (F, /3, t e , 9 e , t max ) dFd(3dt max , (16) 

where e hs(F, /3, t e , 9 e , t max ) is the event detection efficiency (counting efficiency). 

In an actual survey, the optical depth sensitivity, like the efficiency, has to be evaluated 
in a Monte-Carlo simulation. Events are generated, and two histograms in the binned 
(t e , 9 e ) plane are constructed: one counting the events detected by the selection algorithm in 
the given observing conditions, the other including all events above the reference threshold 
A. Weighting each event by F(j)(F)w(j3)S(j3) and summing over events, the two histograms 
evaluate the numerator and denominator of eq. fll6"l) . The sensitivity is then their bin-by-bin 
ratio. 

To illustrate the dependence of the optical depth sensitivity on the luminosity function, 
we consider a simple expression for the detection efficiency. We count all events for which 
the flux increase stays above the threshold A for a duration longer than t\ but shorter than 
£2- In addition, in this example we neglect photometric errors and finite-source effects (so 
to have a dependency on t e only). This translates into e Q b s = 1 for t\ < 2t e ^JP^^ F — f3 2 < t 2 
and e Q bs = otherwise. We calculate the optical depth sensitivity for the luminosity 
functions mentioned previously, and plot the results in fig. 2. Each panel corresponds to 
a value of the threshold magnitude m t h r , and £2/^1 is fixed at 100. Although the details 
differ from luminosity function to luminosity function, the sensitivity curves are similar, 
especially considering that with a more realistic expression for the detection efficiency these 
curves would be smoothed out and that in a comparison with a model they would be folded 
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with the model t e distribution and then integrated. Again, a measurement of the luminosity 
function would reduce even these small systematic uncertainties. 

We finally comment on finite-source effects. With sources of finite angular diameter, 
the flux increase in a microlensing event does not reach the value of the pointlike case 
whenever the lens is over the disk of the source. The exact light curve depends on the 
luminosity distribution on the disk of the star, and in general it is a complicated function. A 
simple account of finite-size effects consists in stopping the flux increase when u(t) < 6 s /26 e , 
where 9 S is the angular radius of the source. Given the minimum detectable flux increase 
A, events with F8(8 s /9 e ) < A but FS(f3) > A are undetected. The efficiency becomes 
a function of the Einstein radius 9 e , as does the sensitivity. The optical depth sensitivity 
as a function of t e and 9 e can be calculated from eq. (IB) once a relation between stellar 



radius and luminosity is specified. We have checked that also in this case the behavior of 
the sensitivity is almost universal. 



5. CONCLUSIONS 

We have proposed an estimator of the optical depth in pixel microlensing experiments 
that involves only measurable quantities. Our new estimator is not limited to events with 
large magnifications at maximum and applies equally well to satellite and ground-based 
observations. 

Our estimator is based on the idea of characterizing the event not by the Einstein 
time and the maximum magnification but by the full event duration at half-maximum and 
the flux increase at maximum. The latter are easily determined from a fit to the event 
light curve. 

The new estimator of the microlensing optical depth is an efficiency-weighted sum over 
events of the products of the fractional flux increase at maximum and of the fraction of 
observing time spent above half-maximum. 

In the ideal case of no noise and complete time coverage, our proposed optical depth 
estimator does not depend on the source luminosity function. In a real situation, an 
experiment is only sensitive to events above a threshold and to a range of Einstein times and 
radii. We have shown that the sensitivity and the effective flux depend on the luminosity 
function in an almost universal way through the surface brightness fluctuation magnitude. 
Residual uncertainties can be removed, if needed, by measuring the luminosity function in 
deep exposures. 
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Our new estimator applies not only to pixel microlensing data but also to blended and 
individually monitored sources. In the latter case, it could be interestingly compared with 
the standard estimator based on Einstein times. If agreement is found, the optical depth 
estimator we propose may indeed open the way to the determination of physical parameters 
from pixel microlensing data. 

The author acknowledges useful comments and discussions with several participants to 
the 4th Microlensing Workshop in Paris, January 1998, where part of this work was first 
presented. 
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Fig. 1. — Dependence of the effective surface magnitude /j, e s on the threshold magnitude m t h r 
for various luminosity functions, fi is the actual surface magnitude and fn is the fluctuation 
magnitude. The lower panel shows the relative systematic error in the effective flux with 
respect to the mean of the luminosity functions considered. 

Fig. 2. — Optical depth sensitivity as a function of Einstein time t e for a simple detection 
efficiency and various luminosity functions. Each panel corresponds to a different threshold 
magnitude m t h r - 
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LUMINOSITY FUNCTIONS 

local [V band] 

47 Tuc [ V band] 

M31 bulge [K band] 

galactic bulge [I band] 




-2 0246 -2 0246 -2 0246 



